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ABSTRACT 

The equivalent width (EW) of the Lya line is directly related to the ratio of star for- 
mation rates determined from Lya flux and UV flux density [SFR(Lya)/SFR(UV)]. 
We use published data -in the literature EW and SFR(Lya)/SFR(UV) are treated as 
independent quantities- to show that the predicted relation holds for the vast majority 
of observed Lya emitting galaxies (LAEs). We show that the relation between EW 
and SFR(Lya) /SFR(UV) applies irrespective of a galaxy's "true" underlying star for- 
mation rate, and that its only source of scatter is the variation in the spectral slope of 
the UV continuum between individual galaxies. The derived relation, when combined 
with the observed EW distribution, implies that the ratio SFR(UV)/SFR(Lya) is de- 
scribed well by a log-normal distribution with a standard deviation of ~ 0.3 — 0.35. 
This result is useful when modelling the statistical properties of LAEs. We further 
discuss why the relation between EW and SFR(Lya) /SFR(UV) may help identifying 
galaxies with unusual stellar populations. 
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1 INTRODUCTION 

Lyman a (hereafter Lya) equivalent width represents a fun- 
damental quantity in the study of Lya emitting galaxies. 
The equivalent width (EW) of the Lya line emitted by 
galaxies is a sensitive indicator of the initial mass func- 
tion (IMF) or gas metallicity from which stars form (e.g. 
■ ISchaerej2 002. 2003). The existence of large equivalent width 
(rest-frame EW > 240 A) Lya emitters has led to specula- 
tion on whether population III f ormation from pristine gas 
may actually have been observed iM alhotra & Rhoads 2002; 

Ijimenez fc Haimanl liooil : iDiikstra fc Wvithel 1200711 How- 
ever, this interpretation depends sensitively on the details 
of radiative tr ansfer through both the intergalactic medium 
(IGM, see e.g. ISchaererll2008l . for a discussion), and the in- 
terstellar medium (ISM, e.g . iNeufeldl Il99ll ; lHansen fc Ohl 

l200fj ; iFinkelstein et al]l2008l ). It is safe to conclude that at 
present, no truly convincing candidates for population III 
galaxies exist. The search for large equivalent width Lya 
emitters - and population III galaxy formation - is a key 
science driver for the observational community, while mod- 
elling Lya radiative transfer in large EW emitters is a chal- 
lenge for theorists. 

The main goal of this paper is to draw attention 
to the fact that Lya EW is frequently discussed inde- 
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pendently from the quantity SFR(Lyq)/SFR(UV) (e.g 
Aiiki et all I2003J ; iFuiita et al.)l2003l; IVenemans et all 12004 



Shimasaku et all 120061 : iGronwall et alj l2007t iTapken et all 
20071 : iPentericci et al.l |2009). This quantity denotes the 
ratio of the star formation rates derived from the ob- 
served Lya flux and rest-frame UV flux density. One can 
consider the quantity SFR(Lya)/SFR(UV) as an "alter- 
native" measurement of EW, because the quantities EW 
and SFR(Lya)/SFR (UV) are directly related (see § I2TT1 of 
this paper, and e.g. IDiikstra fc Wvithel 120071 I Rauch et al 



20081 . iDaval et al.ll200Sl . lNagamine et alJl2008l . lNilsson et al 
2009). Since EW represents a fundamental property of Lya 
emitting galaxies, a more detailed investigation of its rela- 
tion to the ratio SFR(Lya)/SFR(UV) is warranted. 



We derive the relation between EW and 
SFR(Lya)/SFR(UV), and compare with observations 
in § [2] We discuss the implications of our results in § [3j 
Throughout this paper we denote the rest frame equivalent 
width by REW, and the observed equivalent width by 
OEW. The two are related by REW=OEW/(l + z). When 
we write "EW" this refers to both OEW and REW. 
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2 THE CORRELATION BETWEEN EW AND 
SFR(Lya)/SFR(UV) 

2.1 The basis of the Correlation 

The star formation rate derived from the Lya line is ob- 
tained from the better-calibrated Ha star formation indi- 
cator given by SFR(Ha)=7.9 x 10 -42 L Ha Mp yr -1 wher e 
Lho is the Ha luminosity in erg s -1 (e.g. lKennicutt|[l998h . 
For case-B recombination the corresponding Lya luminos- 
ity of the source is Lt, ya = 8.7Ln a , and thus SFR(Lya) = 
9.1 x 10 -43 L LyQ Mq yr" 1 . 

Similarly, the star formation rate obtained from the 
UV luminosity density is g enerally given b y SFR(UV) = 
1.4 x IQ- W L V , UV M Q yr" 1 (|Kennicuttlli99Sh . L„, uv is the 
UV luminosity density in ergs -1 Hz -1 measured at A = 
Auv- To be consistent with published work we shall assume 
Auv = 1400 A. Formally, the st andard conversio n applies 
when A UV = 1500 - 2800 A (e.g. lKennicutt|[l998h . and the 
assumed star formation calibrators may not be completely 
accurate. This only introduces a minor systematic uncer- 
tainty in the connection of UV flux density to an actual star 
formation rate, and our working assumption Auv = 1400 A 
does not affect our main conclusions at all. 

For the star formation calibrators discussed above, the 
ratio SFR(Lya)/SFR(UV) is 



SFR(Lya) 



SFR(UV) 1.4 x 10- 28 L„,£/y 
The Lya rest-frame equivalent width is defined as 



REW 



A L y 

^Ly a L u } Ly a 



(i) 



(2) 



with La, Lya (Li/, Lya) the flux density in ergss -1 A -1 
(ergs s -1 Hz -1 ) of the continuum just redward of the Lya 
emission line, and ALya = 1216 A (^Lya = 2.46 x 10 15 Hz) 
denotes the wavelength (frequency) of the Lya transition. 
We further used the identity XL\ = vL v to obtain the right 
hand side of Eq [2] 

We combine Eq Q] and Eq [2] and find 



SFR(Lya) _ REW L„, Lya 
SFR(UV) ~~ REW C L v ,uv 



(3) 



with REWc 



= 76 A. The ratio of the 



g.lXlO-^LyQ 

continuum flux densities at A = 1216A and Auv = 
1400 A (L^Lyc*/ L v> uv) depends on the slope of the con- 
tinuum. This slope is usually denoted by the parameter 
f3 ee -dlogL A /dlogA (i.e. L\ tx X~ 13 and L„ tx i/ 3-2 ). We 
rewrite Eq|3]as 



SFR(Lya) _ / REW \(v h 



)( 



0-2 
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REW 



(rkw 



(4) 



SFR(UV) VREW c /\i/uv 

Eq g] shows that the ratio SFR(Lya)/SFR(UV) is deter- 
mined uniquely by REW and /3, because REW C , u^ ya , and 
vtjv are constants. The precise choic^] of star formation cal- 
ibrator enters entirely through the value of REW C , and thus 



For example, another conversion factor that is often found 
in the literature is SFR(UV)= 1.25 X IQ~ 2S L V , UV Mq yr -1 
llMadau et al.|[l998l) . Additionally, the conversio n from Lya lu - 
minosity to SFR depends on gas metallicity (sec Schacrcr 2003). 



the slope of the correlation. Provided that the same star for- 
mation calibrators are applied consistently to an ensemble 
of galaxies, the actual star formation rates in these galaxies 
are irrelevant to the existence of the correlation. 

Furthermore, any scatter in the correlation between 
REW and SFR(Lya)/SFR(UV) enters entirely through scat- 
ter in (3. Lya emitting galaxies at z ~ 5.7 show tha t 
< P < 2.4 for > 90 % of the galaxies l|Tapken et alj|2007h . 
This is consistent with the median value /3 me d ~ 1.4 for Lya 
emitti ng objects at z ~ 3 determined by e.g. lVenemans et al.l 
(2005). Thus, we conservatively adopt /3 = 1.2 ± 1.2, which 
translates to G = 0.89to.it for A UV = 1400 A. 



2.2 Processes that affect the Correlation 

Obtaining LLy Q and L Vt uv from the measured quantities 
Fhya (the Lya flux in erg s -1 cm -2 ), and f v ,uv (the UV 
continuum flux density in erg s -1 cm -2 Hz -1 ) is non-trivial 
as it requires various "corrections" that account for the fact 
that only a fraction of all emitted UV and Lya photons are 
observed. None of these corrections are trivial: 

(i) For example, the IGM is expected to be opaque (trans- 
parent) to photons that were emitted blueward (redward) of 
the Lya resonance. To first order, this implies that only 
~ 50% of the Lya photons are transmitted to the ob- 
server. However, peculiar velocities of intergalactic gas (e.g. 
iDiikstra et "all 120071 ; llliev et al.l 120081 ). and radiative trans- 
fer effects in the ISM of galaxies (e.g. lAhn et al.l 120031 ; 
IVerhamme et al . 2008) make this correction quite uncertain. 
The escape fraction of Lya photons from the ISM of galaxies 
is heavily regulated by the pr esence and distribution of dust, 
as we l l as HI gas kinema t ics dKunth et al. 1998; Havcs ~et al.l 
l2009l ; lOstlin et al] |2009| ; lAtek et all 12008 ). This can cause 
the Lya escape fraction to vary significantly between ob- 
jects, or even between different sight-lines within the same 

J II ^1 

object (jLaursen et al . 2009). 

(ii) Dust also affects the UV continuum. The spectral 
slope in the UV conti nuum correlates w i th the amount 
of du st extinction (e.g. ICalzetti et all lT994; Hcc kman et all 
1998). The UV continuum slope for unobscured star form- 
ing galaxies, /3i nt , varies between — 2.6 depending on the 
age, initial mass function, and the st ar formation history of 
a gal axy (see e.g Figs 31 and 32 of iLeitherer fc Heckmanl 
Il995h . Dust lowers the observed UV continuum slope to 
(3 = /3 int - A/3, in which A/3 ~ 1.0[E(B-V)/0.3] (Calzetti 
et al. 2000; where E(B-V) denotes the colour excess). 

Dust reduces the overall U V flux density by a fac tor of 
~exp(-E(B-V)/0.1) (see e.g. IVerhamme et all l200&t ). The 
measured slope of the UV continuum /3 and /or the colour 
excess E(B-V) can constrain the amount of extinction of 
the UV continuum by dust. However these constraints 
are uncertain and depend on the intrinsic UV continuum 
and/or on the precise shape of the extinction curve (see e.g. 
IVerhamme et al.ll2008l ). 



Combined these effects introduce large uncertainties in the 
conversion from observed Lya flux to intrinsic Lya lumi- 
nosity, and from observed UV flux density to intrinsic UV 
luminosity density. Of course, corrections for dust and/or 
the IGM should affect the left and right hand side of Eq|4] 
equally. In practise however, the UV continuum at A ~ AL yQ 
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Figure 1. Comparison of the predicted correlation (Eq[4] solid line) and its 1 cr uncertainties (grey region) with published data (see 
text). Errorbars on all data points denote la uncertainties. Where errorbars are missing, the paper contained insufficient information 
to compute them. We show the relation SFR(Lya)/SFR(UV) = C (REW/REW C ) as the black solid line, where the grey shaded region 
denotes the 1 a uncertainty in C. The grey dashed lines indicate the locations of galaxies that would have spectra with slopes j3 = [—5, 5] 
with steps of 1. The vertical thi ck solid line denotes the maximum REW max = 240 A that can be emitted by 'normal' star forming galaxies 
(e.e. lMalhotra &: Rho ads 2002). Eq[4]converts this maximum REW into a maximum ratio SFR(Lya)/SFR(UV), which is indicated with 
the horizontal thick solid grey line, and its 1 cr uncertainties are indicated with the horizontal thick dotted lines. The majority of data is 
consistent with the correlation given by Eq. [4] Statistically significant outliers correspond to galaxies with an unusual spectral slope of 
their continuum flux density, f3. Right Panel: Same as left panel but focusing on the data points with published uncertainties on both 
REW and SFR(Lya)/SFR(UV). 



is often measured directly from the spectrum, while the con- 
tinuum at A = Atjv is determined from broadband imaging. 
Similarly, the Lya flux that enters the left and right hand 
side of Eq [3] may be obtained differently, e.g. from either 
the spectrum or from narrowband imaging. Clearly, this can 
introduce different systematic uncertainties to the measure- 
ments of SFR(Lya) /SFR(UV) and EW. These uncertainties 
relate to, e.g. slit losses when taking spectra, dependence 
on aperture size and seeing for imag ing, the precise filter 
curve of the narro wband filter (see e.g. lHaves fe Ostlinll2006l ; 
lHaves et al]|2009l . for extended discussions on complications 
that arise when determining f3 and REW from broad and 
narrow band measurements). 



Despite these complications, the quoted uncertainties 
on SFR(Lya)/SFR(UV) and EW should reflect the system- 
atic uncertainties, and we generally expect the data points 
to be consistent with Eq[4] Any statistically significant de- 
viation implies that (i) some corrections are not applied 
consistently between the two measurements [e.g. the IGM 
is corrected for when determining the EW, but not when 
determining SFR(Lya)], (ii) uncertainties in either EW or 
SFR(Lya)/SFR(UV) have been underestimated, or (Hi) the 
continuum of the galaxy has an unusual spectral slope. 

The first two possibilities suggest that Eq 3] provides 
a convenient sanity check to whether the uncertainties in 
SFR(Lya)/SFR(UV) and EW are estimated properly. The 
alternative, an unusual spectral slope is of great scientific 
interest. For example, objects that are dominated by nebular 
emission, such as galaxies that con tain population III stars 
|Schaerer| [2002) or cooling clouds (|Diikstral [20091 ), may be 
dominated by the two-photon continuum at 1216 A < A < 
1600 A (for which (3 < 0). 



2.3 The Correlation in Existing Data 



We show the REW of LAEs versus the ratio of the SFRs 
for several surveys in Figure fTJ The redshifts o f the surveys 
range fr om z ~ 3 — 6.5 . The d ata are from Fuiita et al.l 
J2003l). [Taniguchi et all (|2005l ), IShimasaku et all (1200^)7 
lOuchi et al.l (l2008h (prey da t a poin ts), lAiiki et al.l (|2003l , 
red points), [ Dawson et al.l (|2004l . green points), and 
IVenem ans et al. (2004|, blue p oints). We represent the data 
from IShimasak u et alj |2006) with two sets of points. We 
derive (i) the Lya flux from the narrowband filter, and the 
rest-frame UV flux density from the z-band, and (ii), the 
Lya flux from the spectrum, and the rest-frame UV flux 
density from the i-band. We use either the published values 
of REW and SFR(Lya)/SFR(UV), or - when these are not 
readily available - computed their values from the available 
data. In theory, the correlation in Eq. [4] is independent of 
redshift. However, there are too few measurements available 
to observationally confirm this. The left panel of Figure [1] 
shows that the vast majority of galaxies are consistent with 
the relation given by Eq. [5] and appear to be normal star 
forming galaxies. 



The right panel of Figure [T] highlights the data 
points with published uncertainties in both REW and 
SFR(Lya)/SFR(UV) from the left panel. These points illus- 
trate more clearly that there are no significant outliers from 
the correlation given by Eq[4] We performed a least-squares 
linear fit to these data points of the form log sfr^Jv) = 



a + 6 log 



and obtained a = -0.07 ± 0.22, and b = 



REWe ' 

1.03±0.11. This corresponds to /3 ~ 0.8±3.6. Stronger con- 
straints on /3 are not yet possible given the uncertainties in 
the data. 
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3 DISCUSSION AND CONCLUSIONS 

We investigate the relation between REW and the ratio of 
star formation rates derived from Lya flux, and rest-frame 
UV flux density [SFR(Lya)/SFR(UV)] (EqQJ. This relation 
derives directly from the definition of equivalent width and 
star formation rate conversion factors, and its only source 
of scatter is the variation in the slope of the UV contin- 
uum at ALya < A < Auv between individual galaxies. The 
correlation exists regardless of the assumed star formation 
calibrators (which themselves depend on the assumed IMF 
and gas metallicity) , or the true star formation rates of these 
galaxies. 

Despite their fundamentally tight relation, Lya REW 
and SFR(Lya)/SFR(UV) are often discussed as indepen- 
dent quantities. We investigate their correlation in existing 
data, and find the vast majority of galaxies to be consis- 
tent with the predicted relation (see Figs[T||. The existence 
of the relation has interesting applications, which are dis- 
cussed next. 



3.1 An Empirically Constrained Lya Based Star 
Formation Indicator 

A SFR derived from the UV flux density is likely more reli- 
able than a SFR derived from Lya flux. Lya scatters through 
the ISM and IGM which makes it hard to determine the 
amount of extinction. Our relation can be used to derive 
a more accurate Lya based star formation calibrator, if we 
require that -statistically- the Lya derived SFR should be 
equal to the UV derived SFR. We introduce the constant M 
such that 

SFRx(Lya) =M X SFR(Lya) = SFR(UV). (5) 

The constant M ensures that the SFR derived for a cer- 
tain galaxy from its measured Lya flux is equal to that de- 
rived from its UV flux density. According to Eq U M — 
[C"(REW/REWc)]~ 1 ; this implies that a galaxy with an un- 
usually large REW has M <C 1. Without this correction one 
would overestimate the SFR from the Lya flux alone. 

In the most general case, the probability P(M)dM that 
M lies in the range M±dM /2 is given by (see Appendix lAl|) 

P(M)dM = dMN f P(/3)P(REWm) REW /^^ df3, (6) 

with M the normalization factor which ensures that 
[™P{M)dM = 1, and KEW M (P) = REW C /(MC). Fur- 
thermore, P(REW)dREW denotes the probability that a 
LAE has an observed REW in the range REW±dREW/2; 
P(/3)d/3 denotes the probability that a LAE has an observed 
P in the range P ± dp/2. 

Figure [5] shows the probability distribution P(A4) ob- 
tained from Eq|6] (black solid line). In this calculation we as- 
sume: (i) P(P)dp is a Gaussian with p = 1.2, and 073 = 0.9. 
This c hoice ensures that ~ 90% of the galaxies have < P < 
2.4 (cf. lTapken et al.ll2007h ; (ii) P(REW)dR EW is an expo- 
nenti al with a scale length of REWl=76 A (|Gronwall et al.l 
120071 ). based on observed Lya emitting galaxies at 2 = 3.1 
with REW> 20 A, i.e P(REW)dREWtx exp(-REW/REW L ) 
for REW>REW m in ~ 20 A, and P(REW) = otherwise. In 
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Figure 2. The black solid line shows the probability distribu- 
tion for M =SFR(UV)/SFR(Lya) as given by Eq g] in which 
P(REW) and P(J3) are determined from the observations (see 
text). The red dotted lines mark the 68% confidence interval cen- 
tered on the median M = 1.13 (indicated by the blue dashed 
line). We thus find that(68% of LAEs have M = l.lto*. 
The data points denote the directly measured distribution of 
SFR(UV)/SFR(Lya) collected from the papers used in our anal- 
ysis. The errorbars denote Poisson uncertainties in these data. 
The observed data distribution agrees well with our derived PDF. 
Furthermore, the derived PDF is reproduced quite well by a log- 
normal distribution which is indicated by the grey dot-dashed 
line, which is easily adopted when modeling statistical properties 
of LAEs. 



Appendix lA2l we show that the precise shape of the distribu- 
tion is not very sensitive to the assumed probability density 
functions (PDFs) of P and REW. 

Figure [2] shows that P(M) peaks at M ~ 0.6. The 
function P(M) is highly asymmetric and its median is 
A4 = 1.13, which is indicated as the blue dashed line. That 
is, SFR(Lya) < 0.88 SFR(UV) for 50% of the galaxies (simi- 
larly, SFR(Lya) < SFR(UV) for ~ 57% of the galaxies). The 
red dotted lines mark the 68% confidence interval centered on 
the median value. We therefore find that for 68% of LAEs 
M = 1 .l^o 7 • in other words, the Lya derived SFR lies 
within a factor of ~ 2.5 from the UV derived SFR for 68% 
of LAEs. Overplotted as the data points is the measured dis- 
tribution of SFR(UV)/SFR(Lya) collected from the papers 
that were used in our analysis. The errorbars denote Poisson 
uncertainties. This observed distribution agrees quite well 
with our derived PDF. The grey dot-dashed line indicates a 
log-normal distribution 

r^A^jA., 1 T 1 / P»g.M] - x\ 2 l dM 

P{M)dM = ^m^[-2 V a )\m^W' (7) 

with x = 0.04 and a = 0.35. This log-normal distribution 
provides a decent fit to the derived and observed distribu- 
tion, and is easily adopted when modeling statistical prop- 

erti es of LAEs. 

lOuchi et al.l l)2008h conclude that between z = 3 and 
z = 6 the observed REW PDF is consistent with no redshift 
evolution. This implies that our derived PDF for A4 also ap- 
plies at redshifts greater than z = 3. On the other hand, the 
measured scalelength of the e xponential REW di stribution 
at z = 2.3 is REWl = 48 A (|Nilsson et alj|2009h . causing 
the PDF to broaden and shift to larger M (see Appendix 
IA2[) . Table [1] summarizes the redshift dependence of the pa- 
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Table 1. Fit Parameters for the log- normal PDF for 
P(\ogM)dM. 





X 


a 


z ^3^ 


0.04 


0.35 


^^3~ 


0.20 


0.31 



a We assume that the REW-PDF, and hence the A4-PDF, 
does not to evolve between 2 = 3 — 6 (Ouchi et al. 2008). 
We take the observed 2 = 3 REW-PDF from Gronwall et al. 
(2007) for 2 £3, and the observed z = 2.3 REW-PDF from 
Nilsson et al. (2009) for z <;3. 



probe to cooling clouds or primordial galaxies than the 
spectrum alone. That additional probe is important espe- 
cially at high redshifts, where the IGM may t ransmit only a 
small fraction of the Lya emitted by galaxies (|Diikstra et al.l 
2007). In this case even those star forming galaxies contain- 
ing population III stars may have REW<REW ma x- This 
strongly suggests that the determined REW alone is not 
enough to identify a population III galaxy. 
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rameters describing the log-normal distribution for M. In 
both redshift bins the standard deviation of the log-normal 
distribution is very similar, with a ~ 0.3 — 0.35. The abrupt 
change of the value of x at z — 3 is clearly a crude ap- 
proximation of the real redshift evolution of the 7M-PDF. 
We have chosen this parametrization because it corresponds 
to the simplest description that is consistent with existing 
data. 

In theoretical work the inverse problem often arises 
in which a Lya luminosity must be obtained from a 
physical SFR. In such a case one may write L^ ya = 
l.ixio^crgs- 1 ( M |™ _ t ) and adopt our PDF for M. This 

allows one to assign an empirically calibrated, variable Lya 
luminosity to galaxies of a given SFR. This prescription is 
clearly not perfect, because not all galaxies that are actively 
forming stars show a Lya emission line. For example, only 
20-25% of Lyman Break galaxies (LBGs) at z — 3 have a 
Lya line with REW >,2 0% A and would classify as a LAE 
fe.g. lShaplev et aill2003l ). and the connection between LBGs 
and LAEs is not well understood. Nevertheless, our sug- 
gested prescription for assigning Lya luminosities to galaxies 
of a given SFR is more realistic than the often used one-to- 
one relation L Lyct = 1.1 x 10 42 erg s" 1 SFR/(M Q yr" 1 ). 

3.2 Outliers in the SFR(Lya)/SFR(UV)-REW 
Plane: Signposts for Unusual Galaxies? 

"Normal" star forming galaxies can emit a maximum 
REW m ax = 240 A. Galaxies with REW>REW max may sig- 
nal the p resence of a galaxy that contains population III 
stars (e.g. lMalhotra fe Rhoadsll2002l ) . Using Eq||]this corre- 
sponds to [SFR(Lya)/SFR(UV)] max = 2.8l°; 4 . The major- 
ity of observed galaxies that have quoted uncertainties on 
SFR(Lya)/SFR(UV) are consistent with this upper limit. 

Objects where nebular emission d ominates, such as 
galaxies contai ning population III stars (ISch aerer 20021) or 
cooling clouds <|Diikstrai r2009), may be dominated by the 
two-photon continuum at 1216 A < A < 1600 A. This re- 
sults in unusually negative values for f3. These objects may 
have been identified in the spectrum itself, because devi- 
ations from the correlation are caused by unusual spec- 
tral slopes. However, reliable measurements of the contin- 
uum just redward of the Lya line and at Atjv = 1400 
A provide a long baseline in wavelength, which may more 
clearly reveal the presence of a continuum dominated by 
two-photon emission. This suggests that outliers in the REW 
- SFR(Lya)/SFR(UV) plane may provide a more sensitive 
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APPENDIX A: THE PROBABILITY 
DISTRIBUTION FUNCTION (PDF) FOR M 

Al Derivation of P(M)dM 

We employ the notation of probability theory, where the 
function p(y\b) denotes the conditional probability density 
function (PDF) of y given b. The PDF for y is then given 
by p(y) = J p{y\b)p(b)db, where p(b) denotes the PDF for 
b. Similarly, we can write p(y) = J J p(y\b,a)p(b)p(a)dbda, 
where p(y\b, a) denotes the conditional probability density 
function (PDF) of y given 6 and a, and where we assume that 
a and b are independent. We can thus write the probability 
P(M)dM that M lies in the range M ± dM/2 as 



P{M)dM = dM M I d(3 dREW 

J~oo ./EW mi „ 

P(M\(3, REW)P(REW)P(/3), 



(Al) 



with REW m in the minimum REW, and P(M\l3, REW) the 
conditional PDF M giveifl f3 and REW. 

Because we know that J °° P(M\f3, REW) AM = 1, and 
that for a given (3 and REW there is only one solution A4 = 
[C(REW/REW C )] _1 (Eq [5)1, we write P(A4|/3,REW) = 
5 U (M - [^(REW/REWc)]" 1 ) = <5 D (ff[REW]), where 5 D (x) 
denotes the standard Dirac-delta function. We eliminate the 
integral over REW by using the property of Dirac delta func- 
tions that 8n(h[x]) = j$t^t, where x Q denotes the (real) 



root of h{x). If we denote the root of g(REW) with REW.M 
then |s'(REWai)| = RE w and obtain EqU 
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Figure Al. This Figure shows the dependence of P(_M)dM on 
various model parameters that describe the PDFs for (3 and REW, 
which include j3, op, REW m i„ and REWl- Each curve represent 
a calculation in which one of the model parameters - as indicated 
in the top right corner of the Figure - was varied. 
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Figure A2. Same as Fig lAll but for model parameters indicated 
in the top right corner. 



A2 Dependence of P(M)dM on Model Parameters 

We investigate the dependence of P(M)dM on various 
model parameters that describe the PDFs for (5 and REW, 
which include /3, 073, REW m i n and REWl- Figures IAT1 and 
IA2I show the PDFs that we obtain when varying one of the 
model parameters as indicated in the top right corner of each 
Figure. For example, the red dashed line shows the PDF that 
we obtain for (3 = 0.0. 

Figures IA1I and IA2I show that the PDF does not vary 
significantly when changing the parameters ap, REW m in- 
The PDF appears to be most sensitive to the scale length 
of exponential REW distribution. The value REWl = 48 
A correspond s to t he scale length that was derived by 
iNilsson et ail (2009) for 2 = 2.3 LAEs. At these lower red- 
shifts there is a lower fraction of large EW emitters, which 
pushes the A4 PDF to larger values. 



2 We assume that REW and /3 are independent variables. In 
the case that larger data sets demonstrate that this assump- 
tion is false, then one needs to replace P(REW)P(/3) with 
P(REW|/3)P(/3) in EqlATl 



